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ABSTRACT 

We propose a new cosmological test of gravity, by using the observed mass fraction of X-ray 
emitting gas in massive galaxy clusters. The cluster gas fraction, believed to be a fair sample 
of the average baryon fraction in the Universe, is a well-understood observable, which has 
previously mainly been used to constrain background cosmology. In some modified gravity 
models, such as f{R) gravity, gas temperature in a massive cluster is determined by the ef¬ 
fective mass of that cluster, which can be larger than its true mass. On the other hand. X-ray 
luminosity is determined by the true gas density, which in both modified gravity and ACDM 
models depends mainly on Ob/Um and hence the true total cluster mass. As a result, the stan¬ 
dard practice of combining gas temperatures and X-ray surface brightnesses of clusters to infer 
their gas fractions can, in modified gravity models, lead to a larger - in f{R) gravity this can 
be 1 /3 larger - value of Ub/Om than that inferred from other observations such as the CMB. 
A quick calculation shows that the Hu-Sawicki n = 1 f{R) model wit h | f rq I = 3 5 x 10~^ 

is in tension with the gas fraction data of the 42 clusters analysed bv lAllen et al.l (l2008h . We 
also discuss the implications for other modified gravity models. 
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1 INTRODUCTION 


In recent years, atte mpts to understand the origin of the accel erated 
cosmic expansion teiess et alJll998l : IPerlmutte r et al.lll 99^ hav e 
led to a large number of theoretical models Icopeland et al.l200^ . 
Apart from the current standard model, which assumes that the ac¬ 
celeration is caused by a cosmological constant A (hence the name 
A-cold-dark-matter, or ACDM), these model can be roughly put in 
two categories: dark energy - which replaces A by some dynamical 
field, and modified gravity - which assumes that there is no exotic 
matter species beyond the standard CDM model but gravity is not 
desc ribed by General Relativity (GR) on cosmological scales (see, 
e.g.. lClifton et al.l2012l : ljovce et al.l2015l . for some recent reviews). 
Some of the models in the l atter class, such as the chameleon theory 
dKhourv & Weltmanl2004l) . of which the well-known f{R) gravity 
dCarroll et al.l2005ti is an example, have been active research topics 
in recent years. 

Ultimately, any new cosmological model or theory of gravity 
should be put to rigorous tests against observational data. For this 
reason, a number of tests have been proposed or applied in the past 
to examine the viability of the models (see, e. g.. lUombrisej[20l^ : 
IPe Martino et al.ll2015f . for some recent reviews in f{R) gravity). 
The present paper shall follow the same line to propose a test using 
observations of galaxy clusters. 

Galaxy clusters are the largest gravitationally bound and viri- 
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alised objects in our Universe. The most massive clusters observed 
today typically have masses in the range of ~ Mq, 

of which the dominant component is dark matter. These are homes 
to galaxies, stars and eventually lives, which together hold the vast 
majority of the information that can be extracted from cosmological 
and astrophysical observations. In dark energy or modified gravity 
theories, the different cosmic expansion histories and gravitational 
laws between part icles can have sizeable effects on h ow the clusters 
form and evolve. ISchmidt. Vikhlinin & Hid d2009l) : ICataneo et al.l 
d2014h . based on this observation, have placed constraints on f(R) 
gravity using cluster abundance data. In theories such as f{R) grav¬ 
ity, massive and massless particles feel different strengths of grav¬ 
ity, thus allowing these theories to be constrained b y comparing the 
so-called dynamical and te nsing masses of clusters dSchmidi2010l : 
ZhaOj^i_^Koyai^ 201 ih. Combin i ng wi th lensing observations, 
Terukina et al.l d2014 ): IWilcox et al.l d2015h obtained even stronger 
constraints on f{R) gravity. 

Given the abundant information associated with these rich ob¬ 
jects, one expects that they will provide a wealth of other potential 
tests of new cosmological models. The arrival of the era of preci¬ 
sion cosmology lends this perspective both more interest and more 
support. In this paper, we propose to utilise the observationally in¬ 
ferred mass fraction of hot X-ray emitting gas in galaxy clusters as 
a new test of gravity. 

The gas fraction in clusters, /gas, is a well established and un¬ 
derstood observable in the standard cosmological model, which can 
be used to place strong constraints on background cosmology (see. 
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e.g.. lAllen et al.ll^04l2008l . and following-up works, for some ex¬ 
amples). The basic assumptions (or approximations) are (i) clusters 
are in hydrostatic equilibrium between thermal pressure and grav¬ 
ity and (ii) as the largest objects in the Universe, the cluster baryon 
fraction, dominantly contributed by gas, is a faithful representation 
of the cosmological average baryon fraction in which fib 

and flm ar e respectively the f r actional mass density of bary ons and 
all matter iWhite et al.fl 19931 : lEke, Navarro & Frenklll99^ . Using 
(i), one can find, from the observed X-ray temperature and surface 
brightness profiles, the mass profiles of baryons and all matter in¬ 
side a cluster, and consequently /gas (r) - the profile of gas fraction. 
Combined with (ii), this can have a say about flm provided that fib 
is measured elsewhere (e.g., from big bang nucleosynthesis or the 
cosmic microwave background (CMB)). 

In modified gravity theories, the hydrostatic equilibrium inside 
clusters is changed as a result of the different law of gravity. Hence, 
a cluster can have a higher dynamical mass, with the same baryonic 
mass inside, leading to a that is lower than the true /gas which 
is related to fib /flm. If a cosmologist wishes to infer fib /flm from 
fgta in a modified gravity universe, a correction has to be done to 
the end result, which can lead to inconsistency with other observa¬ 
tional determinations of fib /flm, and therefore a constraint on the 
gravity theory. We shall demonstrate the potential constraints from 
this test using f(R) gravity as an example. 

This paper is organised as following; in Sect.|^we briefly de¬ 
scribe the f(R) gravity theory and its equations which will be used 
in the discussion below. In Sect. [3] we give a more detailed account 
of the physics related to the gas fraction test described above. Then 
we present a numerical example in Sect.|4l which shows how cur¬ 
rent data of cluster gas fraction alone can give powerful constraints 
on gravity. We discuss the results and their implications in Sect.j^ 


2 THE F(R) GRAVITY THEORY 


This section is devoted to a quick overview of f(R) gravity. It will 
be kept brief and only include essential equations, given that there 
are already many p apers in the literatu re covering this topic. 

f(R) gravity dCarroll et al.l2005h is a simple generalisation of 
GR, by replacing the Ricci scalar, R, in the Einstein-Hilbert action 
with an algebraic function f(R) 


S = j [R + f{R)]+Cr, 


( 1 ) 


where Mpi is the reduced Planck mass, = SttG with G being 

Newton’s constant, g is the determinant of the metric and Cm 
is the Lagrangian density for (normal plus dark) matter fields. The 
model is defined by specifying the functional form of f{R)- 
The action in Eq. Q leads to a modified Einstein equation 


G^1/ f rR^u 




- V^V./k = SttGT” , (2) 


in which G^i^, R),,,, are respectively the Einstein and Ricci tensors, 
fn = d//di?, the covariant derivative compatible to the metric 
□ = V“ Va and is the matter energy momentum tensor. 
Eq. I2I1 can be considered as the standard Einstein equation of GR 
with an extra scalar field, fn, whose dynamics is governed by 


n/ii ^ 3 “ /rR + 2 / + 87rGpm), 


(3) 


where pm is the mass density of baryons and dark matter. As we are 
interested in late times, photons and neutrinos will be neglected. 


On scales well inside the Hubble radius, and for the models to 
be considered, it is safe to wo rk with the quasi-static approximation 
( iBose. Hellwing & Lij2015l) . in which the scalar equation becomes 

V^/fl = [Rifn) - R + SttG [pm - Pm)] , (4) 

where V denotes the three dimensional gradient, a is the scale fac¬ 
tor, and an overbar takes the background value of a quantity. Notice 
that R can be expressed as a function of fn by inverting /^(i?). 
Similarly, the modified Poisson equation in this limit reads as 

V"(E> = (pm - Pm) + ^a^[R (/fl) - i?] , (5) 

where is the Newtonian potential. 

Eq.© implies two limits of the behaviour of f{R) gravity: 

(i) When /_r is small, or more accurately, when |/_r| <C |4?|, it 
recovers the well-known GR solution R = Syr G/Jm, and so Eq. © 
reduc es to GR as well. This is the chameleon jKhourv & WeltmM 
l2004h regime which any viable f{R) model must be in to pass the 
stringent solar system and terrestrial tests of gravity. 

(ii) When | /_r | ~ G (| <& |), the second term on the rhs of Eq. © 
is negligible compared with the first term, so that we have a gravity 
that is 1/3 stronger than in GR. This is ususally known as the non¬ 
chameleon, or unscreened, regime. 

It is evident that the unscreened regime mostly happens where 
is shallow, or in extensive regions of low density. On large scales, 
matter density is close to the cosmological average, and so the total 
gravity in f[R) gravity is enhanced within scales comparable to the 
Compton wavelength of the scalar field fu (which in most models 
of interest is in the range G(1 ~ 10)/i“^Mpc). This naturally leads 
to an enhanced large-scale structure formation, and features such as 
over-abundant and more massive galaxy clusters - a topic which has 
been extensively studied previously. This will also be the topic that 
we focus on in this paper. 


3 CLUSTER GAS FRACTION 

Galaxy clusters are the largest bound objects in the Universe, whose 
masses are dominated by the dark matter component, with the bary¬ 
onic masses dominated by X-ray emitting intracluster gas, which is 
heated to temperatures of the order of keV during virialisation. It is 
the mass fraction of this gas component that we will employ to test 
the theory of gravity here. 

In this section, we shall first give a brief overview of how the 
baryon fraction can be estimated observationally, and how it can be 
used to constrain cosmological models and their parameters. Then, 
we will discuss how this process might be affected if the underlying 
theory of gravity is modified. For simplicity, we shall neglect other 
baryonic components than the intracluster gas in our analysis unless 
otherwise stated. 


3.1 The standard ACDM model 


In the standard cosmologic al scenario, halo density profiles can be 
universally described by the lNavarro. Frenk & Whita ( 1 19971 . NEW) 
fitting formula, which is often expressed as 


p{r) 


_Ps_ 

(r/rg) (1 + r/ra)^ ’ 


( 6 ) 


in which p{r) is the halo mass density as a function of the distance, 
r, from the halo’s centre, pa is a characteristic density and Vs is the 
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scale radius. We shall assume that the halo is spherically symmetric 
and well relaxed throughout the analysis, unless otherwise stated. 

The mass of the halo can he obtained by integrating the NFW 
profile from r = 0 to r = Ra , in which Ra is the edge of the halo 
and is defined as the radius within which the average mass density 
is A X pcv\t{z), with pcTit{z) = 3H{zf /8nG the critical density 
at the redshift 2 when the halo is identified. This leads to 


Mhaio = 47rpsr® 


In (1 + Ra/ts) 


Ra/tb 
1 + Ra/^s 


(7) 


Observationally, the total and baryonic masses of a cluster can 
be obtained by measuring its X-ray surface brightness profile, and 
the temperature profile of its X-ray gas. For a dynamically relaxed 
system that consists of dark matter and baryonic gas, a hydrostatic 
equilibrium can be achieved, which satisfies the following equation 


1 d 
pgas(r) dr 


Pgas (r) 


GMtot(< r) 
^2 


( 8 ) 


in which Mtot{< r) is the total mass of dark matter and gas within 
radius r, and Pgas(r), Pgas{r) are respectively the density and pres¬ 
sure of the gas at r. For simplicity, we neglect non-thermal pressure 
in our discussion (the effects of non-thermal pressure, however, are 
taken into account in the error budget when modelling the relation 
between /gas and cf. Eq. ( I19t below). 

For an ideal thermal gas, its pressure and density are related to 
its temperature, Tgas, as 


prup 


-PgasTgi 


(9) 


where k is the Boltzmann constant, and the mass and number den¬ 
sities of the gas particles are connected by pgas = ^mptigas, where 
trip is the proton mass and p the mean molecular weight. Applying 
Eqs. {Hig, we obtain 


GM{< r) fcTgas(r) 

r prup 

which, if evaluated at the halo edge (r = Ra), gives 

GMhalo _ kTgss{RA) 

Ra pmp 

In the mean time, the X-ray emission of galaxy clusters is pro¬ 
duced mainly by thermal bremsstrahlung radiation, leading to a X- 
ray surface brightness profile that also depends on the gas density 
and temperature profiles, pgas(r) and Tgss{r). In terms of gas den¬ 
sity, the gas mass can be expressed as 

rr 

Mg^s{<r)= / 47rr'^drVgas(r-0 oc Pgas(0)ro, (12) 

^0 

where Pgas(O) is the gas density at the cluster centre, and we have 
assumed that Pgas(r) = Pgas(0)p(r/ro) with g{x) some function 
describing the profile and ro a characteristic scale. In the isothermal 
/3 model, for example, we have 


d In pgai 
dlnr 


+ ■ 


d In Te, 


dlnr 


( 11 ) 


dlnpgas(r) ^ d In Tgas (r) 


dlnr 


dlnr 


( 10 ) 


Pgas(r) 


Pgas(O) 



(13) 


where is a dimensionless constant. _ 

The bolometric luminosity is given by (e.g.. lSasa Ol 19961) 

Lx{<r)cx[ 47rr'^dr'rg\4^(r')pgas(r') oc Pgas(0)ro, (14) 
Jo 

where we have neglected the proportionality coefficients which are 
irrelevant here. 


Combining Eqs. GlIIl, we have 
Mgas oc oc oc (1 -f 


(15) 


where dn and dx are respectively the luminosity and angular diam¬ 
eter distances and are related by dn = (1 -F zfdx- Here we have 
used Lx = 4:TTd^fx where fx is the X-ray flux, and ro = dxQo, 
in which 0o denotes the angle spanned by ro at redshift z. Note that 
fx (or equivalently the surface brightness) and 0o are the observed 
quantities in this description. In real observations, one has the sur¬ 
face brightness and temperature profiles, or equivalently Lx{< r) 
and Tgss{r), using which Eqs. HOI 1141 can be solved simultane¬ 
ously to find Mhaio(< r) and Mgas(< r) (note that the innermost 
regions of clusters are often excluded due to co mplicated processes 
such as cooling flow. For example, the study of lAllen et al.l ( 1200 Sh 
uses only the data within 0.7 ~ 1.2i?2500 to measure /gas). These 
provide the necessary information to find 


/gas(r) 


Afgas(< r) 
Mhalo(< r) ■ 


(16) 


Often in observations, people quote the value of /gas at r = R 2500 
as the cluster gas fraction. We note in passing that here Mhaio enters 
the picture only through its gravitational effect on gas particles, as 
this fact is important for the discussion of modified gravity below. 

To see how this can be used to constrain background cosmol¬ 
ogy, we note that, once the right-hand side of Eq. (ED is known by 
observations, we also have a fixed numerical value of the left-hand 
side and hence have the relation 


Mhaio OC Ra OC dx- 


(17) 


Eqs. (ED ED imply that the measured cluster gas fraction depends 
on redshift and angular diameter distance in the following specific 
way: 


/gas(z) = oc (1 -F zfd!f^. 


(18) 


As mentioned earlier, we expect the cluster gas fraction to be a rea¬ 
sonably fair sample of the mean cosmological baryon fraction and 
is therefore roughly independent of 2 for massive clusters at low 2 . 
This is a reasonable assumption which should hold regardless of the 
cosmological model/parameters. The measurement of the apparent 
gas fraction, as described above, involves the angular diameter dis¬ 
tance dx and is indeed dependent on both the cosmological model 
and its parameters. Therefore, with incorrect cosmological models 
or parameters, the constancy of the true gas fraction (hereafter /gas) 
is not guaranteed to be re flected in the obse rved value fg^ ( 2 ). This 
provides a powerful test dAllen et al.l2004h of background cosmol¬ 
ogy and can be used to constrain cosmological parameters. 

Were galaxy clusters perfectly fair samples of the average mat¬ 
ter components in the Universe, their baryon fraction would just be 
flb/Hm- The true situation, h owever, is more co mplicated. To take 
the complexiti es into accou nt, Allen et aD(l2008h . improving on the 
earlier work of lAllen et alj ( l2004h . propose the following model of 
the relation between cluster gas fractions and Hb/flm, 


K-yhiz) Hb 
1 -F s{z) ’ 


(19) 


in which: 

• K is a. constant accounting for systematic effects such as the 
calibration of inst rument and X-ray m odelling - it is assumed to be 
K=l.G± 0.1 in lAllenetalJ ( l2008h : 

• 7 models the non-thermal pressure support in galaxy clusters 
which can cause a bias in the estimate of fg^ of about 9%; 
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• b{z) = ho{l-\-a\,z) is the so-called ’depletion’ factor which 
is inspired by the observa tion that the baryon fraction at -R 2500 in 
non-radiative simulations llEke. Navarro & Frenklll99ill is actually 
smaller than fib/f^m, with 60 = 0.83 ± 0.04 and an small indicat¬ 
ing a weak redshift evolution below z = 1\ 

• s{z) = so(l+crsa) accounts for the fact that a small fraction 
of baryons can be in the form of stars, with sq = 0.16 ± 0.048 and 
—0.2 < Qa < 0.2 describing its redshift evolution. 

The model in Eq. ( 119b indeed has a weak redshift dependency. 
However, any additional dependency from the observed /gas would 
imply that one is using the wrong background cosmology to extract 
data, cf. Eq. ( fTSb . This is, to be clear, in the framework of standard 
GR. 


3.2 Modified gravity scenarios 


In many modified gravity theories, including f{R) gravity, the way 
in which the trajectories of massive test bodies - e.g., galaxies, stars 
and gas particles - respond to the underlying matter distribution is 
different. This change of the dynamics of test bodies is sometimes 
described as the change of the dynamical mass of matter. Massless 
particles, such as photons, coul d behave differently : in some theo¬ 
ries, s uch as the Galileon model ^Nicolis_etaj200^ De^vet et al.l 
1 20091) and the K-mouflage model ISrax & Valageas.l2014al f9). nho- 
tons can also feel a different ma ss of matter, but in other models, for 
example f{R) gravity and the lPvali, Gabadadze & Porratil l l200Ct 
DGP) model, photon trajectories depend on the matter distribution 
in essentially the same way as in GR, since the conformal coupling 
does not affect geodesics of massless particles. For distinction, the 
mass felt by photons is usually called the tensing mass. The differ¬ 
ences in the dynamical and lensing masses of galaxy clusters hav e 
been used to constrain fiR) gravity in, e.g.. lTerukina et al.l ( 120141 ): 
I Wilcox et all 120 ISt) 

As mentioned above, gas particles, like dark matter particles 
and galaxies, do feel the dynamical mass of a cluster. In f{R) grav¬ 
ity, the same cluster can have a dynamical mass 4/3 times its value 
in GR. This maximum enhance factor of 1/3, however, is not nec- 
essarily realised in all clust ers, because of the chameleon screening 
iKhourv & Weltm^l2004l) . The screening helps to reduce the dif¬ 
ference between the dynamical and lensing masses, especially for 
more massive clusters. Consequently, constraints relying on the dy¬ 
namical masses of clusters are in general weaker than those coming 
from astrophysical considerations. Nevertheless, they have cleaner 
physics than that of astrophysical observables - which can often de¬ 
pend on whether the considered astrophysical system lives inside a 
screened cluster - and ar e amongst the tightest constraints obtained 
using cosmological data iTerukina et al.ll20M) . 

Evidently, without a reliable measurement of the lensing mass, 
it is difficult to tell whether one observes a cluster of true mass M 
in GR, or one with a smaller mass in f{R) gra vity, since both ha ve 
the same dynamical masses. This has motivated iHe et al] ( l2015ah to 
propose the concept of effective haloes. Briefly speaking, the idea 
is to redefine the right-hand side of the modified Poisson equation 
in f{R) gravity, Eq. (|5}, so that it can be rewritten as 


V^<I> = 47rGo^(5pm,eff. 


( 20 ) 


In the above, <1? is the Newtonian potential that determines the dy¬ 
namics of massive bodies and Spm = Pm — pm is the density pertur¬ 
bation of (dark plus baryonic) matter. 5pm,eff is the effective density 
field, with which the Poisson equation takes exactly the same form 
as in GR [cf. Eq. (I20b1. In this way, the complicated new physics in 
f{R) gravity is absorbed into pm.eft, and with that solved (e.g., in 


numerical simulations) one c an in principle pr oceed assuming GR 
as the true theory of gravity. iHe & Lil ( l2015bh . for example, show 
with hydrodynamical simulations that cluster gas temperatures de¬ 
pend only on the masses of the corresponding effective haloes, and 
that with certain rescaling depending on the effective halo mass the 
scaling relations - such as the Lx-M relation with Lx the X-ray 
luminosity - in f{R) gravity can be derived reliably using existing 
knowledge of GR. 

In the left panel of Fig.[T] we present the gas temperature pro¬ 
files for standard haloes in GR and effective haloes in f{R) gravity, 
both in the mass bin 10^® ~ Though there are dif¬ 

ferences in the inner regions - which could be due to different halo 
density profiles or scree ning - we notice t hat beyond ~ 100h“^kpc 
the two agree very well. iHe & Lil ( l2015bh find that the average gas 
temperatures in the two also show very good agreement, and indeed 
the temperature-mass relation is barely distinguishable in the two 
models, provided that effective haloes are used in f{R) gravity. 

The fact that the cluster gas temperatures depend on the mass 
of the effective haloes is as expected, since for relaxed systems the 
virial temperature depends on the Newtonian potential, which does 
not distinguish between standard (GR) and effective (,f{R) gravity) 
density fields. For a polytropic gas with an equation of state Pgas oc 
Pgas in which the constant F ^ 1, the hydrostatic equation implies 
that (see, e.g.. lMo. van den Bosch & Whitell2()ll1) the temperature 
can be analytically expressed as a function of the potential 4?. 

Let us now consider two haloes, one identified in the standard 
dark matter field in a model with GR as the gravity theory, another 
from the effective density field in f{R) gravity. The profiles of the 
two haloes are the same so that Eq. ( 120b sees no difference in them. 
Since gravity only enters the picture through Eq. dlOb , we make the 
following two observations/predictions: 

• the gas temperature profiles are the same in these two haloes; 

• the two sides of the spherical hydrostatic equation, Eq. Oil, 
are the same for the two haloes. 

Since Tgas (r) and d In Tgas (r ) /d In r in Eq. dlOb are the same, 
we conclude that d In pgas (r) /d In r is also the same in the haloes. 
This, however, does not necessarily mean that the haloes have iden¬ 
tical gas density profiles, because we can rescale pgas by a constant 
factor without changing d In pgas (r)/d In r. To confirm this, in the 
middle panel of Fig.[T]we compare the gas density profile of f{R) 
effective haloes with that of ACDM haloes of the same mass, and 
find that the two show a constant shift by 1/4 beyond r ~ 
kpc. Note that the simulated haloes in the plot do not have perfectly 
identical total - standard or effective - mass profiles, which is why 
in the middle panel of Fig.[T]the red crosses and red dashed line do 
not agree on scales below ~ 100/i“^kpc. However, as mentioned 
above Eq. ( 122b . in real observations, such innermost regions are not 
used in the determination of /gas anyway. 

As a result, to obtain the gas density profile, we need further, 
independent, information to fix its normalisation - as opposed to 
its shape - which brings us back to the measurements of cluster X- 
ray surface brightness. Inspecting the equation for the cluster X-ray 
luminosity, Eq. {H, we notice that the luminosity density (i.e., the 
integrand) depends on (i) the physical gas density Pgas, and (ii) the 
gas temperature Tgas which, as we have seen above, depends on the 
total mass of the effective halo. Consequently, should the physical 
gas densities be the same for f{R) effective and ACDM haloes of 
the same mass, there would be no difference in their X-ray surface 
brightness profiles. 

However, despite the standard (GR) and effective (f{R) grav¬ 
ity) haloes above having identical gas temperature and halo mass, 
their actual (physical, or lensing) masses are different, and it is im- 
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Figure 1. (Colour online) Left panel, the halo gas temperature profiles from the non-radiative hydrodynamical simulations of iHe & lH i2Q15bl) for the standard 
ACDM model (black solid line) and f{R) gravity with |/ho| = 10“^ (blue dotted line for standard haloes and red dashed line for effective haloes. The same 
line styles ai‘e used for the other two panels). Middle panel: The same as the left panel, but for the gas density profiles. The red crosses are a rescaling of the 
ACDM curve by 3/4 to take into account the fact that for a ACDM halo and a f{R) effective halo of the same mass, their true masses differ by 1/3. Right 
panel: The halo gas fraction profiles from the same simulations; the red crosses are the ACDM result scaled by 3/4, which is almost identical to the red dashed 
line, cf. Eq. E). The profiles shown here are stacked results of the haloes in the mass range ~ h ^Mq, for illustration purpose; we’ve checked 

that results from other halo mass bins follow the same trends. 


portant to remember that characterises the amount of the actual 
mass in the Universe. If, as we have assumed so far, the gas fraction 
in clusters is a fair sample of the cosmological value, it would be the 
ratio of the gas mass and actual halo’s mass that satisfies Eq. d- 
Gas fractions inferred ohservationally, in the way described in the 
previous subsection, are in fact the ratio of the gas mass and that of 
the effective halo. If we denote the ratio of the effective and actual 
masses of a halo by rj, then 1 ^ ^ 4/3, depending on the actual 

mass and environment of the halo, its redshift, as well as the f{R) 
model parameterl3- Here, as we are interested in the most massive 
clusters, with Mhaio ~ Mq, we can for simplic¬ 

ity neglect the impacts of the halo’s environment, so that rj mainly 
depends on z, i.e., rj = ri{z), for a given halo mass. 

Consider the extreme case in which 77 = 4/3 as example. The 
apparent /gas inferred from X-ray cluster observations would be 


robs _ .^gas 

^ Mhalc.efl 


3 Mgas _ 3 

4 flThalo,actual 4 


( 21 ) 


This is also confirmed by hydro simulations, as shown in the right 
panel of Fig.[T] In that plot, we see that if haloes are defined using 
their actual masses in f{R) gravity, then they share the same /gas as 
ACDM haloes of the same mass (black solid vs. blue dotted lines). 
For effective f(R) haloes, on the other hand, their /gas profiles are 
a constant downward shift by 1 /4 from the results of ACDM haloes 
of the same masses (red dashed line vs. red crosses), which is just 
what Eq. (ED predicts. 

It is worthwhile to pause for a moment and try to understand 
the physics behind the behaviour of Fig. [D It may seem surprising 
that, although the gas density profiles in ACDM and f{R) standard 
haloes are significantly different within r ~ 100 /i“^kpc, their gas 


^ More accurately speaking, effective haloes are identified from the effec¬ 
tive density field, Pm,eff in Eq. (20), and standard haloes are identified from 
the physical density field pm ■ They do not necessarily share the same phys¬ 
ical particles. Here, for simplicity, when talking about the effective and ac¬ 
tual masses of some halo, we mean the masses of the effective and standard 
haloes that would be considered as matched haloes in the two catalogues. 


fraction profiles are very close to each other. This suggests that the 
two also have different dark matter (or total) mass profiles to cancel 
the differences in pgas(r). In other words, pgas(f) follows pDM(f), 
which is the physical dark matter mass density, in the same way in 
ACDM and f{R) standard haloes, even though they have different 
mass and even more different potential profiles. To understand this, 
we note that in Eq. (d the G on the left-hand side and Tgas on the 
right-hand side are both modified in f{R) standard haloes. Indeed, 
if one assumes hydrostatic equilibri um and Pgas oc the n the 
gas density profile can be written as (Komatsu & SeliaMl200lh : 

Pga;^(f) _ r - 1 G/r trip Mhaio C f^m{x) 

pr-i(o) - r i?vufcTgas(0) m(c) 7o *2 ^ ’ 

in which m(r/ra) = M(< r)/4'Kpsrf Eq. ( 122b is derived under 
the assumption of self-similarity of gas density profiles, but the key 
point therein, that the modified gravity effects on G and Tgas can be 
cancelled out, is not affected by this assumption. If this cancellation 
happens, then the gas density profile is determined by the total mass 
profile under the assunmtion of hydrostatic equilibrium, regardless 
of the theory of gravitjQ 

In clusters, gas is heated by accretion shocks during the assem¬ 
bly of the halo, a process which involves the conversions of energy 
from gravitational to kinetic (that of the cold accreted gas) and then 
to thermal (via shocks). Assuming a complete thermalisation, the 
post-shock gas temperature is proport ional to v?„r„u, with Uinfaii the 
infall speed of the accreted gas (e.g., van den Bosch & Whitel 
l201lh . Consequently, energy conservation implies that the final gas 
temperatures in the central regions will be affected in the same way 
as wfnfaii of the cold gas and hence G in modified gravity. Of course, 
this is only an approximation, and the cancellation of the effects of 
modified gravity on G and Tgas depend on various factors includ¬ 
ing the screening and formation history of a cluster, which is not ex¬ 
pected to be complete. However, Fig.[Dsuggests that it works pretty 


^ Note in Eq. (22) it is pgas(^)/pgas(0) that is determined by m(r/rs). 
The normalisation of pgas will then be fixed by the total gas fraction inside 
the halo. 
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Figure 2. The inferred /Om (triangles with 1 a error bars) using the for the 42 clusters from lAllen et alj ilQQSh . as a function of the cluster redshift. 

4 cases are shown, with different underlying models of gravity: GR (upper left), and Hu-Sawicki n = 1 f(R) model with |/ro| = 3 x 10“® (upper right), 

5 X 10~® (lower left), 10“^ (lower right). The horizontal solid and dotted lines are respectively the mean and Icr range of fJb/rim from Planck CMB data. 

A good match between /gas and CMB for the GR case, and progressively worse matches for the f{R) models, can he seen by a quick inspection by eye. For 
simplicity, all haloes are assumed to have a mass of 7.5 X and a concentration of 3.3 when determining the effect of the chameleon screening. 


well for the haloes we use here. We checked explicitly that it works 
slightly less well for more massive haloes, for which the agreement 
between the /gas(f) in ACDM and f{R) standard haloes is slightly 
less perfect - this may be because those haloes became unscreened 
only very recently. 

The argument above in theory also applies to effective haloes, 
for which G is the same as in GR, but the effects of modified gravity 
are incorporated in Mhaio- However, in the effective halo case the 
normalisation is different because of the different total gas fraction 
(see footnote 2) - although the shape is the same - hence the nearly 
constant rescaling of the red dashed curve compared with the black 
solid one in Fig.[T] 

Coming back to the discussion prior to the previous three para¬ 
graphs, our result suggests two possible tests of f{R) gravity: 

(i) If an observer actually lives in a universe shaped by f{R) 
gravity, then the true cluster gas fraction is given by ~ Vfg^s- 
Assuming that Eq. Gl) still holds for /gas, the observer will need 


to do the following transformation to get the true S7b /Hm: 


robs 1 K'yb{z) 

Cb 


Cb 



fib 

7) 1 + s( 2 ) 


true 

12 m 

true 

= V 

12m 


As (Hb/Gm)oba depends only on the actual observational data, the 
observer will obtain the same value as an observer in a standard GR 
universe would do. The resulting (Gb/Gin)true might then be too 
large to be compatible with other constraints, such as the one from 
the CMB. 

(ii) Alternatively, if one takes the S7b /Cm measured by other 
probes as the true value and starts from there, then Eq. | |2H implies 
that the observed cluster gas fraction f°^ will be smaller than what 
the ACDM model and simulations predict. Because of the time de¬ 
pendence of r]{z) (see above), if the f{R) model parameters hap¬ 
pen to take the values for r) to evolve from 1 to 4/3 between z = 1 
and the present for the clusters of interest, there may also be an ap¬ 
parent decrease of /g^/ ( 2 ) as 2 decreases, by a maximum of 25%. 
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Table 1. The assumed ranges and distributions of the model p arameters in 
Eg. jl 91. For more details the readers are referred to Sect. l3.1l or lAllen et alj 


param 

physical effect described 

mean ± stddev 

prior 

K 

overall calibration 

1.000 ± 0.100 

Gaussian 

7 

non-thermal pressure 

1.050 ±0.050 

Unifonn 

bo 

gas bias: normalisation 

0.825 ±0.175 

Uniform 

oib 

gas bias: evolution 

0.000 ±0.100 

Uniform 

^0 

stellai* fraction: nonnalisation 

0.160 ±0.048 

Gaussian 

Ots 

stellai* fraction: evolution 

0.000 ± 0.200 

Unifonn 


4 NUMERICAL EXAMPLES 


In this section, we use a simplified example to illustrate the power 
of the cluster gas fraction test proposed above. F or this, we will use 
the gas fraction data of the 42 clusters studied bv l Allen et alJjlOO^ 
Table 3). As described above, these /gas data are obtained by fitting 
the gas temperature and X-ray surface brightness profiles of these 
clusters simultaneously, assuming NFW profiles for the total mass 
in clusters. We have also found that, in the context of f{R) gravity, 
as long as we use effective haloes, the dynamics of gas particles can 
be calculated using standard gr avity theory. There fore, in this work 
we can directly take the data o/ lAllen et alJ l l2008tl as f°^i, bearing 
in mind that the cluster mass inferred therein would be the effective 
mass and therefore can be different from for unscreened 
clusters, cf. Eq. HD. 

To obtain an estimation of the mean and standard deviation of 
(flb/f2m)o|3s = K7i,(z) fgas from each cluster, random samples of 
size 10® are drawn for each parameter or data: K, 7 , bo, an, so, cts 
a nd Of these, is taken, for a given cluster, from Table 3 
of lAllen et alJ HoOSti . and is assumed to satisfy a Gaussian distribu - 
tion with mean and standard deviation given bv l Allen et al.l HoOSli . 
The other parameters and their distributions are shown in Table [D 
We therefore obtain 10® realisations of from which 

its mean and standard deviation can be calculated. This procedure 
is repeated for all 42 clusters. 

The estimation of the effects of modified gravity, i.e., the fac¬ 
tor 77 ( 2 ), is more complicated, since it depends on the cluster mass, 
density profile, environment, redshift, as well as the f{R) parame¬ 
ters. Because the main purpose of this paper is to illustrate the basic 
idea, we shall leave a full analysis using real cluster data for future 
work, and instead adopt a simplified modelling. The cluster masses 
are assumed to be the same, with Mhaio = 7.5 x Mq, for 

all 42 clusters, since this is a typical value for massive X-ray clus¬ 
ters. The cluster’s radius (R 200 or 7?vir) is taken to be 1.5/i“^Mpc 
and its concentration parameter, c = i? 20 o/fs (or Rvii/rs) is 3.3. 
The cluster is assumed to live on the cosmological background, so 
that the ratio of its effective and actual m asses can be approximated 
as (see, e.g.. lLi. Zhao & Kovamall^l2h 


fz) = 


-^VThalo.eff 

Mhaio ,actual 


min 



fnC) 

2<E-n ’ 3 J 


(24) 


where fniz) is the background value of fa at redshift 2 , and 


TrAfhalo ln(l -t- c) 

Rhs.io ln(l + c) - c/(l -I- c) 


-5 X 10"®, 


(25) 


is the Newtonian p otential at the edge of an NFW halo. We adopt 
the f{R) model bv iHu & Sawickil 120071) with n = 1, for which 


fnif 


dm “t“ 4dA 
dm(l + ■ 2 )® -f 4dA 


fm- 


(26) 


The value s dm = 0.316 and Da = 0.684 are tak en from the latest 
results of IPlanck Collaboration. Ade et all ItOISIi . As a result, the 
physics of modified gravity is completely governed by /ho, which 
is the present-day value of fa. Once this is specified, we can obtain 
77 ( 2 ), and therefore infer (db/dmltr^g given (db/dmloba * of 
a cluster. 

In Fig.lDwe show the (db/dm)tjue result obtained from 
for 4 different cases: standard ACDM (upper left), and f{R) grav¬ 
ity models with |/ho| = 3 x 10“® (upper right), 5 x 10“® (lower 
left) and 10"'^ (lower right). For comparison, we have also, in each 
panel, plotted the mean value (solid line) and Icr confidence lev el 
(dotted) of db/dm from IPlanck Collaboration. Ade et aP HolSiFl 
The results from the 42 clusters, with Icr errors, are shown as sym¬ 
bols. 

A quick naked-eye inspection shows that the /gas method and 
the CMB observation give compatible db/dm if one assumes the 
ACDM paradigm (upper left). The f{R) model with |/ho| = 10"'* 
(lower right), on the other hand, leads to a significantly higher value 
of db /dm than what CMB says, and is therefore inconsistent. The 
other two cases are more interesting: for |/ho| = 5 x 10"® (lower 
left), 77 ( 2 ) increases to 4/3 at 2 ~ 0.45, while for | /ho | = 3x10"® 
(upper right) 77 ( 2 ) only increases to ~ 1.28 at 2 = 0. In both cases, 
however, the inferred values of db/dm are still substantially larger 
than the Planck result, especially for the low -2 clusters. This shows 
that cluster gas fraction can be a potentially powerful test of gravity, 
using X-ray observations only. In such tests, lensing data can be a 
useful addition, but is not necessary. 

The test can be done in an alternative way. For this, we assume 
the value of db/dm obtained from CMB observations, and check 
what value of an observer would have found if living in a /(7?) 
universe. The idea is that, if this v alue differs too m uch from what 
our observers have told us (e.g., in lAllen et al.l200^ . then it would 
place a constraint on the extent to which the assumed f{R) model 
can deviate from standard ACDM. As in the previous case, we draw 
random samples of size 10® for the parameters K, 7 , 60 , an, so, cts 
from which we find 10® realisations of Kyh{z)/{1 -\- s{z)). Then, 
by modelling modified gravity effects using Eq. HD, we compute 
the mean /gas for the 4 models shown in Fig.Hl and these are shown 
as curves in Fig. together with the observed values of /gas from 
lAllen et aP HOO^ Again, we note that current data favour ACDM 
over all three variants of f{R) gravity. As clusters are less screened 
at late times, we find that low -2 data is more useful in constraining 
the model than high -2 data. 


5 DISCUSSION AND CONCLUSIONS 

In this paper, we proposed a new cosmological test of gravity, by 
inferring the cosmic baryon fraction from the apparent gas fractions 
of massive clusters, and comparing with the results from other, less 
model-dependent, measurements such as the CMB. In theories with 
a stronger gravity, the apparent gas fraction is smaller than that in 
ACDM for a fixed Db/Dm. Reversely, if the observed value f°^ is 
fixed, we would find a higher value of Db /Dm than in GR, that can 
be inconsistent with the model-independent measurements. Taking 
the Hu-Sawicki f{R) model as an example: our quick calculation 


® Note that the f{R) models studied here have practically identical CMB 
power spectra as the ACDM model with the same Dm and Da . As a result, 
using CMB data only, the constraints on cosmological parameters such as 
Dm would be the same in all these models. Because of this, the CMB con¬ 
straints are less model-dependent. 
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Figure 3. (Colour online) The evolution of /gas in 4 models w i th the same 
cosmic value of fin/f2m from IPlanck Collaboration. Ade et alj jlOlSl) . The 
models are, respectively, ACDM (black solid line), and Hu-Sawicki n = 1 
/(-R) model with |/flol = 3 X 10“® (blue dotted line), 5 X 10“® (green 
dashed line) and 10“^ (red dot-dashed line). Black triangles wi th en'or bars 
are the /gas values of the 42 clusters used in lAllen et al j i20o3. Tab 3). The 
blue shaded region denotes the standard deviation around the mean /gas for 
the model with = 3 X 10“®, for illustration, which shows that the 

theoretical uncertainty is roughly of the same order as cun'ent observational 
errors in /gas; therefore, the constraining power can be further improved if 
either of these uncertainties is reduced in the future. 


shows that model parameters | fjio | ~ 3 — 5 x 10~^ are in tensio n 
with the gas fraction data of the 42 clusters from lAllen et alJ l l2008h . 
though a more rigorous constraint will be left for futur e work. 

/ gas has been a rather widly-used observable (e.g.. lWhite et al.l 
Il993l) . and its power in constraining cosmology - in particular dark 
energy mod els - is conv i ncingly demonstrated in various previous 
works (e.g.. lSasakill996l : rAllen et alj2004l2008h . The inclusion of 
baryons opens a new dimension for tests of gravity, since ultimately 
most cosmological observables can be tracked back to lights emit¬ 
ted by interactions involving baryons. In the mean time, the physics 
of the X-ray-emitting hot gas in massive clusters is relatively clean, 
making it easier both for the modelling and to use the observational 
data. As an example, the assumption that gas temperature depends 
on the gravitational potential and our main conclusions that (i) 

- the true gas fraction - is unchanged with modified gravity while 
(ii) /gas® is changed a re supported by hydrodynamical simulations 
in f{R) gravity (e.g.. lHe & Lill2015bl) . Some uncertainties remain 
in relating /gas to fib /flm, but these have been included in the error 
budget estimate above. Furthermore, within our current state of un¬ 
dersta ndi ng, slightly changi ng its modelling (e.g., from lAllen et al.l 
( I 2 OO 4 I) to lAllen et al.l l l2008h ) does not change results drastically. 

Here, we w ould like to emphasise the use of effective haloes 
jHe et alj2015ah in our analysis. Though the idea has a simi lar ori¬ 
gin as that of the dynamical mass of halo (e.g., ISchmidll2010h . there 
are fundamental differences. Dynamical mass is a certain attribute 
of a given halo which is defined in the standard way, while effective 
halo is a completely new way to define and identify haloes. Given 


an effective halo, all gravitational effect can be calculated from GR, 
and in particular this means that the way in which /gas is currently 
extracted from observational data - and the resulting /gas results - 
can be directly used for our purpose. Thus, with a little extra effort 
from the people who generate a halo catalogue, the analyses of end 
users can be made much more straightforward, and this provides an 
efficient bridge between simulators, theorists and observers. 


One may naturally wonder about the generality of this method. 
As a cosmological test, it relies on galaxy clusters being totally or 
partially unscreened. Because we are talking about massive clusters 
which tend to be better screened, this test, like most other cosmo¬ 
logical ones, will probably not be able to constrain fao to substan¬ 
tially smaller than the quoted values here. However, it does provide 
a fairly clean test - with good observational data available - that has 
the potential to place one of the strongest constraints from cosmol¬ 
ogy on f[R) gravity. Also, o ne can always combin e /gas and other 
observables, such as lensing 1 TemkinaetalfcOl^, cluster scaling 
relations (e.g. [Arnold. Puc hwein & SDringelll2014l) . and cluster gas 
pressure profiles ( De Martino et alj2014h . to place joint, and likely 
stronger, constraints. In principle, the test would be more powerful 
if observational data for smaller galaxy clusters (e.g., those in the 
mass range 10^® ~ Mq) and galaxy groups are included, 

because these objects are less screened and so gravity deviates more 
from GR in general. This, however, requires a better understanding 
of the feedbacks in different models, which are not well studied so 
far. 


The test can be applied not only to f{R) gravity and the more 
gener al chameleon the ory, but also to sim ilar models such as dila- 
tons jBrax et al.ll201lh and symmetrons jHinterbichler & KhourvI 
I 2 OIOI) . These models are all featured by a universal coupling of all 
matter species to a scalar field, that effectively enhances the gravity 
for all particles (at least in unscreened regimes). There are models 
In which only certain matter species, e.g.. dark matter, experiences 
the scalar coupling: therein, baryons can still feel a different gravity 
depending on how the dark matter particle mass evolves with time, 
in which case the proposed test does apply. In addition to these the¬ 
ories, we have mentioned above the DGP, Galileon and k-mouflage 
models. In the first two classes, the de viations from GR are strongl y 
suppressed inside dark matter haloes tearreira et alj|201^l2014a) . 
and so we do not expect the new tes t to work. For K-mouflage, as is 
for the so-called non-local gravity dMaggiore & Mancarellall2()14l : 
IPirian et al.l2014l : lBarreira et alj20i4^ . there can be a time evolu¬ 
tion of Newton’s constant inside clusters, making it possible to use 
this test. However, one needs to bear in mind that in many of these 
theories the background evolution history is also modified, and that 
can affect the /gas test (whether it leads to degeneracies or stronger 
constraints can only be told by a case-by-case study in the future). 

As mentioned earlier, the aim of this paper is to illustrate the 
main idea of using /gas as a test of gravity theories, and therefore 
we have made a simplified estimate and have not quoted any numer¬ 
ical results on the confidence levels of the constrained fao. A more 
complete and rigorous analysis will require one to relax the simpli¬ 
fication that all observed clusters share the same mass, radius and 
concentration, and use the real observational results of these for all 
clusters. If the cluster mass is obtained from its dynamical effects, 
we also need to account for the fact that different clusters may have 
experienced different degrees of screening, and so a more accurate 
modelling of the screening is needed to compute the cluster mass 
profile. These will be left for future work. We note that hydrody¬ 
namical simu lations for modified gravity theories started to appear 
recently (e.g.. lAmold. Puchwein & SDringej2014l : lHammami et al.l 
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l2015l:lHe&Lil2015bh . and such works will be useful for improving 
the constraining power of this test in the future. 
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